An extended class of minimax generalized 
Bayes estimators of regression coefficients 



We derive minimax generalized Bayes estimators of regression co- 
efficients in the general linear model with spherically symmetric errors 
under invariant quadratic loss for the case of unknown scale. The class 
of estimators generalizes the class considered in Maruyama and Straw- 
derman (2005) to include non-monotone shrinkage functions. 
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1 Introduction 

In this paper we consider minimax generalized Bayes estimators of the re- 
gression coefficients in the general linear model with homogeneous spheri- 
cally symmetric errors. We start with the familiar linear regression model 
Y = A/3 + e where Y is an N x 1 vector of observations, A is the known 
N x p design matrix of rank p, f3 is the p x 1 vector of unknown regression 
coefficients, and e is an iV x 1 vector of experimental errors. We assume e 
has a spherically symmetric distribution with a density a~ N f(e'e/a 2 ), where 
a is an unknown scale parameter and /(•) is a nonnegative function on the 
nonnegative real line, which satisfies 
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The problem is to estimate (3. The least squares estimator of (3 is (3 — 
(A'A)~ 1 A'y. In order to treat the estimation problem from the decision- 
theoretic point of view, we measure the loss in estimating (3 by b with so 
called scale invariant "predictive loss" functions 

L(b,f3,a 2 ) =a- 2 (b-P)'A'A(b-P). (1) 

Then the risk function of an estimator b is given by R(b, (3, a 2 ) = E[L(b, (3, a 2 )]. 
The least squares estimator (3 is minimax with constant risk. Therefore, b 
is a minimax estimator of (3 if and only if R(b, (3, a 2 ) < R(/3, (3, a 2 ) for all 
(3 and <r 2 , the search for estimators better than (3 is a search for minimax 
estimators. It is, of course, well known that estimators dominating (3 exist 
when p > 3. In this paper, we study generalized Bayes minimax estimators. 

To simplify expressions and to make matters a bit clearer it is helpful to 
rotate the problem via the following transformation, so that the covariance 
matrix of (3, a 2 (A' A)^ 1 becomes diagonal. Let P be the p x p orthogonal 
matrix of eigenvectors of (A'A)^ 1 , with d\ > d,2 > • • • > d p as eigenvalues, it 
follows that 

P'{A'A)- 1 P = D, P'P = I p 

where D = diag(rf 1? . . . , d p ). Let also Q be an N x N orthogonal matrix such 
that 

' D -l/2p' 





QA 



Next, define two random vectors X = (Xi, • • • , X p )' and Z = (Zi, ■ ■ ■ , Z n )' 
where n = N — p by 

A ' QY. 



Then (X', Z')' has the joint density given by 

<T- p - n f({(x-9y(x-0)+z'z}/<j 2 ), (2) 

where 9 = P'(3. Notice also that X and Z'Z can be expressed as D~^ 2 F "$ 
and (y — A(3)'(y — Af3) 1 respectively. Also, as is customary, denote Z'Z by 
S. We assume throughout this paper that p > 3 and n > 1. 

The original problem is thus equivalent to estimation of 9 under the loss 
function L(5, 9, a 2 ) = (5 — 9)' (5 — 9)j a 2 . We will consider the problem in this 
equivalent canonical form. 
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This paper is best viewed as a companion paper to Maruyama and Straw- 
derman (2005). In that paper a class of minimax generalized Bayes estimators 
was derived for the above canonical problem when the errors were normally 
distributed, and a subclass was shown to be generalized Bayes and minimax 
for the entire class of spherically symmetric error distributions. 

In this paper we enlarge the class of generalized Bayes minimax estimators 
for each of above classes of distributions. To do so we must first enlarge 
the class of minimax estimators to accommodate non-monotone shrinkage 
functions. Section 2 is denoted to this generalization. Section 3 develops an 
extended class of generalized Bayes minimax estimators for spherical normal 
error distributions. To a large degree, the results therein represent extentions 
to the case of unknown variance, of the results of Maruyama (1998, 2004) in 
the known variance case. Section 4 extends these results to the case of general 
spherically symmetric error distributions with unknown scale. A subclass of 
the estimators in Section 3 is shown to be generalized Bayes and minimax 
for the entire class of spherical error distributions simultaneously (subject to 
finiteness of moments). 



2 Minimaxity 

In this section, we give a sufficient condition for minimaxity in the general 
spherically symmetric case. 

Theorem 2.1. Suppose (X',Z')' has a distribution given by (T5]). Then 5^ 
given by 




is minimax if 

^ Un + 2U(w) - 2(p - 2)} - 40'H{1 + <j>(w)} < 0. (4) 
w 

The interesting point of the theorem is that the sufficient condition for 
minimaxity does not depend on /. Such distributional robustness has already 
been noted in the literature, but only the following tractable subset of the 
above result is typically used: 
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Corollary 2.1. If <fi is monotone nondecreasing and < <p(w) < 2(p — 
2)/(n + 2), 5^ is minimax in the general spherically symmetric case. 

For completeness we give the proof. We use the version derived in Kubokawa 
and Srivastava (2001) but earlier versions appear in Robert (1994) and else- 
where. 



Proof. Let 

and define 

E f [h(X,Z)\ 

E F [h(X,Z)} 



F(x) 





h(x,z)a- N f 
h(x,z)a~ N F 



f(t)dt 
[x - 6)'{x - 

cr 2 (J z 

(x-ey(x-e) z 'z\ , , 

^ J -i '- + — dxdz, 



z'z 



H 7T dxdz 



a- 



where h(x, z) is an integrable function. Note that F(x) is not necessarily a 
probability density, and hence E F is not necessarily an expectation symbol 
in any strict sense. 

The identities corresponding to the Stein and chi-square identities for the 
normal distribution (Stein (1973) and Efron and Morris (1976)), 



E'[(Xi - di)h(X, Z)\ = ^[(d/dXiMX, Z% 
Ef[Sg{S)]=a 2 E F [ng{S) + 2Sg'{S)l 



(5) 
(6) 



where S = Z'Z, are useful in our proof. 
The risk of 9^ is given by 



R(9,a 2 



E f 

R(6,a 2 ,X) + E f 



ey^-e)^ 2 
s 2 



o- 2 \\x\ 



\x\ 



s 



2E f 



E 



s 



a- 



X 



-Xi(Xi — 9i 



\X\ 



Let W = \\X\\ 2 /S. For the second term in 



E f 



S_ 

a^\\X\ 



S^ 2 



\x\\ 



E L 



(n + 2) 



using ([oj, 

<f) 2 (W) 



S 

we have 



(7) 



W 
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For the third term in (17j) , using (jSJ), we have 



IXI 
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p- 



W 



+ 2<- 



Hence we have 



R{9,a 2 ,e^) -R(9, a 2 , X) 



E 1 



<ttW) 

w 



{(n + 2)<f>{W) - 2(p - 2)} - 4<j)'(W)(l + <j>{W)) 



which completes the proof. 



□ 



An immediate corollary, which we believe is often more tractable, is given 
by dividing the left hand side of the inequality (jl]) by {1 + 4>(w)}(f)(w)/w. 



Corollary 2.2. (i) 5^ given by (j3J) is minimax if 

(n + 2 )<p(w) - 2(p - 2) _ < Q 



1 + 0(W) 



(8) 



(ii) Suppose < <j>{w) < Mi and w<j)'(w)/<j)(w) > — M 2 . T7ie estimator 5^ 
given by ([3]) zs minimax if 



[n + 2)M 1 - 2(p-2) 
1 + Mi 



+ 4M 2 < 0. 



It is clear that part (ii) of the corollary above allows be non-monotonic. 
Actually we show in the next section that there exists a class of minimax 
generalized Bayes estimators with non-monotone (/>. 



3 An extended class of Generalized Bayes min- 
imax estimators in the normal case 

In this section, we extend a class of generalized Bayes minimax estimators 
in Maruyama and Strawderman (2005) in the normal case. 
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Suppose the sampling distribution of (X', Z')' is normal with covariance 
matrix a 2 I n and mean vector (8', 0')'. As in Maruyama and Strawderman 
(2005), the class of hierarchical priors we consider is as follows; 

%, A ~ N p (0, v -\l - \)\~ 1 I P ), r] ~ V e , A ~ A a (l - X) b . (9) 

Then the marginal density of X, S, A and rj is proportional to 



exp 

rp 



2 \ 11 11 1 — A 11 11 J 2 J (1 - A)p/ 2 " 6 



oc exp (-y(l + Aw)) TjpW+e\p/*+*(i - \)\ (10) 

where w = ||x|| 2 /s. Under quadratic loss, the generalized Bayes estimator 
for such a hierarchical priors can be expressed as (TTOl) . 

= E(r,0\X, S) ( E(\ n \X,S) \ ( 4><W)\ 

e = E(Ksr = i 1 - e(„\x,s) ) x = { 1 -—) x - 

When p/2 + n/2 + e + 2 > 0, 

r/P /2 + n/ 2+e+ i exp + Aw) j rfr? K (1 + Aw) - p / 2 _ n/M; (12) 

hence 

, x f 1 A p/2+a+1 (l - A) 6 (l + wX)~ p/2 ~ n/2 ~ e ~ 2 dt , 

4>(w) = w-^-, , (13) 

/„ Af/ 2+a (l - A) b (l + wA)-p/ 2 -"/ 2 - e - 2 ^ 

which is well-defined for a > —p/2 — 1 and 6 > — 1. 

Here is a result of Maruyama and Strawderman (2005). 

Theorem 3.1 (Maruyama and Strawderman (2005)). Suppose b > 0, e > 
—p/2 — n/2 — 1 and —p/2 — 1 < a < n/2 + e. Then 4>(w) is monotone 
increasing and approaches (p/2 + a + 1)/ (n/2 + e — a) as w — ► oo. Hence 8$ 
with <f> given by f[T3"j) zs minimax if 

p/2 + a + l p-2 
< ^ < 2- 



n/2 + e — a n + 2 
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The key assumption of the theorem is b > while the estimator with 
4>(w) itself is well-defined when b > — 1. In fact, there is a big difference 
between the two cases: b > and — 1 < b < 0. An immediately apparent 
difference is that (1 — X) b is unbounded for — 1 < b < while it is bounded for 
b > 0. Technically, because of this unboundedness, the standard integration 
by parts technique fails to work well when — 1 < b < 0. Furthermore the 
shrinkage factor <b(w) for —1 < b < is not often monotonic, which means 
the tractable sufficient condition for minimaxity given by Corollary 12 .11 is not 
applicable. 

We will see below that proving minimaxity with b < requires a different 
approach from that in Maruyama and Strawderman (2005). Actually it will 
be done through the expression of cb by hypergeometric functions 

F{a ,p-, r ,,) = i + ±^ x 4 (14) 

where (a)i = a • ■ ■ (a + i — 1). The following lemma summarizes the rela- 
tionships we use. All formulas are from Abramowitz and Stegun (1964) and 
the number following AS in each expression below is the formula number in 
Abramowitz and Stegun (1964). 



F(a, (3- 7 ; z) = wm ^, 7) m C - ^"^(l - tz)- a dt (15) 



Lemma 3.1. • AS. 15. 3.1 

n 

r(/?)r( 7 -/3) 7 

when 7 > j3 > 0. 
AS. 15.3.4 

F(a, (3; 7; z) = (1 - z)- a F(a, 7 -#7; z/(z - 1)). (16) 
AS. 15.2.25 

7(1 - z)F(a, (3; 7; z) - 7 F(a, /3 - 1; 7; z) (17) 
+ z(j - a)F(a, P; 7 + 1; z) = 0. 



AS. 15.2. 18 



(7 - a - /3)F(a, /5; 7; z) - (7 - a)F(a - 1, # 7; 2) (18) 
+ /3(l-z)F(a,/3 + l; 7 ;z) = 0. 
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• AS. 15. 1.20 

lim F(a 3- r z) = 1(7)1X7 - - - P) (19) 

when 7 7^ 0, —1, —2, . . . and 7 — a — 3 > 0. 

• AS. 15.3. 10 

lim{- log(l - z)}- 1 ^, /?; a + /?; z) = i^^Jl ■ (20) 
z-*i 1 

• AS. 15.3.3 

F(a, 3- 7; 2) = (1 - zy-^F^ - a, 7 - 8; 7; 2). (21) 

• AS. 15.2.1 

(d/dz)F(a, 3- 7; 2) = {a/3/ 7 }F(a + 1, /3 + 1; 7 + 1; z). (22) 

We first use Lemma [3. II to re-express <p(w). 

Lemma 3.2. Lei be given by f|T3|) . T/ien provided a > —p/2 — 1, b > — 1 
and p/2 + n/2 + e + 2 > ; </>(«;) is expressed as 



<f>(w) = \J. (23) 



1 - G(u 

l+e—a 
p/2+a+l 

where 

F{b, p/2 + n/2 + e + 1; p/2 + a + b + 2; u 



G(«) 



F(6 + 1, p/2 + n/2 + e + 1; p/2 + a + 6 + 2; v) 

Proof. Using (|T5j) and (TT6T) . we have 
1 

^-1^ _ ^7-/3-1 (! + (24) 

= r(/3 r(7)~ /3) (1 + 7 " ^ ; 7; + 2)) 

when 7 > /3 > and hence by ({TBI 



1; p/2 + o + lF(6+l,c;d+l;^) 
0W = 1^ d F(6+l,c;rf;,) ' (25) 
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where v = w/(w + 1), c = p/2 + n/2 + e + 2, and d = p/2 + a + b + 2. Using 
(HZD and ([IB]), we have 

F(6 + l,c-l;d;v) 

0(WJ = -1 + 



-1 



1 - u)F(6+ l,c; d; u) 
c- 1 

n/2 + e - a + (p/2 + a + l)G(u) 
1 - 



26 

□ 

The next lemma gives properties of ratios of hypergeometric functions 
such as G, which we employ in demonstrating minimaxity. 

Lemma 3.3. Let 

KM = Ji^7" ] i (^) 
F(6 + 1,^7; u) 

for (3 >0 and 0<6 + l<7<6 + l + /3. Then, 
(i) K(0) = 1 
(wj lim^i K{v) = 0. 

fwj If b > or if —1 < b < and 7 > /3 ; i^(n) zs monotone decreasing. 

(iv) If — 1 < 6 < and 7 < /3 ; ine minimum of K(v) is a negative 
value and K(v) approaches /rom the below as v approaches 1. yMso 
inf <„<iA» > 6/(6 + 1). 

Proo/. First note that K(0) = 1 by (I27jl . Next by (fT9l and (|2T|), 

limd - ft 7; .) = r(7) r r '°^r 7) < 28 > 

2-1 1 {a)l [pj 



when 7 > and 7 - a - /3 < 0. Hence we have by (j25|i . ( 1201 and ( |T9i) 
respectively, 

■ ( l_^)7-^ r(7W- 7 ) if7 _ 6 _^ <0 

F(6, /3; 7; v) ~ <| {- log(l - ^)}r^& if 7 - 6 - /3 = (29) 

' if0< 7 -6-/3<l 



and also by 



. 6 _ 1H9 r( 7 )r(& + i + /3-7) 



F{b + 1, ft r , v) « (1 - v f-^-P W r( V 6 + 1)r( ^ ~ ^ 7 - & - P < 1 

(30) 

where f(v) ~ g(u) means lim^i f(v)/g(v) = 1. We easily see that the 
ratio F(6, /3; 7; v)/F{b + l,/3; 7 ;u) = i^(t>) goes to zero as t> — >■ 1 under the 
assumption 7 — b — (3 < 1. Hence part (ii) follows. 

When b > 0, K(t> ) is decreasing from the monotone likelihood ratio prop- 
erty of the kernel of k(v). Hence the first assertion of part (iii) follows. 

When — 1 < b < and 7 — f3 > 0, the numerator F(b, (3;^;v) is always 
positive because it can be rewritten as (1 — t>) 7 ~ fc_/3 F(7 — b, 7 — /3; 7; v). Also 
the numerator F(b, /3;7;f) with — 1 < 6 < is decreasing in t> and the 
positive denominator F(b + l,/3;7;t>) is increasing in t>. Hence the second 
assertion of part (iii) follows. 

To show part (iv), note that T(x) > if x > and T(y) < if — 1 < y < 0. 
By assumption -1 < 6 < 0, 7 < (3 > 0, 0<6+l<7<6+l + /5. It 
then follows using the additional assumptions in the first and third lines 
of that there is exactly 1 negative factor in each constant term and 

hence each constant term is negative. Since, also, the denominator of K(v), 
F{b + l,/3;7;t>), is positive, it follows that K{y) approaches from below 
as v approaches 1. Let K{y) take on its minimum value at Vq. Using the 
formula 



(A/B)' = (B'/B){A'/B' - A/B} 
and fl22|) . we have K{vq) = Ki(vo) where 

b F(b+l,8 + l;i + l;v) 
1{V> 6 + lF(6 + 2,/3 + l;7 + l;u)' 

Since 6/(6 + 1) < and 6 + 1 > 0, Ki(v) is increasing in v. Therefore 

K(v) > K(v ) = K^Vq) > K^O) = 6/(6 + 1). (31) 
Hence part (iv) follows. □ 

The following corollary gives the behavior of (j){w) and follows immedi- 
ately from Lemma 13.21 and 13.31 
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Corollary 3.1. Assume that e > —p/2 — n/2 — 1, — p/2 — 1 < a < n/2 + e 
and b > — 1. T/ien 



(%) lim^^oo </>(«;) = a where a = (p/2 + a + l)/(n/2 + e — a), 
fzz) When 6 > min(n/2 + e — a — 1, 0), monotone increasing. 

(Hi) When —l<b< min(n/2 + e — a — 1, 0), <f>(w) is not monotonic. 
(iv) When 

n/2 + e — a 



< b < min(n/2 + e - a- 1,0), (32) 



< <£M < 



p/2 + n/2 + e + 1 

1-6/(6+1) (p/2 + a + r 



l/a + 6/(6 + 1) n/2 + e - a + 6(p/2 + n/2 + e + 1) ' 



Note: By analogy with the known variance case (Alam (1973) and Maruyama 
(2003, 2004)), it may be expected that the choice — 1 < 6 < would lead 
to a non-monotone 0. However part (ii) of the corollary shows that this 
need not be true in the unknown variance case, and in fact monotonicity 
depends on the relationship between n, a and e. The addition of the restric- 
tion fl32l) in part (iv) is necessitated by the fact that lim^-i 6/ (6+1) = 
— oo. Hence a value of 6 close to —1 would cause the upper bound on <f), 
{1 — 6/(6 + l)}/{l/a + 6/(6 + 1)} to be negative. Thanks to the restriction, 
the denominator is positive. 

The next result gives a lower bound for w<p'(w)/<p(w). 

Lemma 3.4. 

£H b/2 + a + 2)6 

4>{w) 6 + 1 K ' 

provided that e > —p/2 — n/2 — 1, —p/2 — 1 < a < n/2 + e — a and —1 < 
b < min(ra/2 + e — a — 1, 0). 

Proof. Note 

{wA(w)/B(w)}'/{A(w)/B(w)} = 1 + wA'(w)/A(w) - wB'(w)/B(w). 
Hence 

w<f/(w) , f /'A^+^a-A^l+wAr^MA 
[ J 1 Ap/ 2 + a + 1 (l-A) b (l + wA)- c rfA 
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where c = p/2 + n/2 + e + 2. Using f|T5|) and f|T6|) . we have 
w4>'(w) _ / | + a + 2 F(6 + 1, c + 1; d + 2; u) 



[ d+1 F(6+l,c;d+l;u 
| + a + 1F(6+ l,c+ l;d+ l;v 



d F(b + 1, c; d; v) " ' ' 

where f — w/(w+ 1) and d = p/2 + a + 6 + 2. Using ffTTl) and ffT8l) . 

^ . (F(b+ l,c+ l;d+ l;u) F(6 + 1, c + 1; d; u) . 



</>(«;) 1 l,c;d+ l;u) F(6 + 1, c; 

= (g + a + 2) ^ c ' d + 1 ' u ) -(g + a + l). 



2 + l,c;d + l;v) v 2 y F(6 + 1, c; d; u) ' 

Note F(b,c;d + l;v) > F(b,c;d;v) since — 1 < b < 0. For t> which satisfies 
F(b, c; d; u) > 0, we have 

- Ftp, c; d; v) ( + * - ^±^-\ > 0, 



<f>(w) " ' \F(b+ l,c;d+ l;v) F(b + 1, c; d;v) 

because F{b + 1, c; d; i>) > + 1, c;d + l;v). For w which satisfies F(b, c; d + 
l;t>) > > F(b,c;d;v), W(p r (w)/<p(w) is clearly nonnegative. For v which 
satisfies F(b, c; d + 1; t> ) < and — 1 < b < min(n/2 + e — a — 1, 0), we have 

**W> ( g + a + 2) f(^+M >(g + Q + 2): & 



~ v 2 l,c;d+ l;v) ~ v 2 '6+1 

by Lemma 13.31 (iv). □ 

The main result is the following. 
Theorem 3.2. Suppose e > —p/2 — n/2 — 1. 
(i) [monotone <fi] When 

- P /2-l<a< c(p '" )( " /2+ ; ) - p/2 - 1 (34) 
1 ~ l + c(p,n) 

where c(p, n) = 2(p — 2)/(n + 2) and 6 > min(ra/2 + e — a — 1, 0), i/ie 
generalized Bayes estimator is minimax. 
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(ii) [non-monotone <p] When 

- P /2-l<a< c( ^" )(n/2+ / e) -/ /2 - 1 (35) 

1 + c(p, n) 

-1 < -(w + 2) , c ^ n )~ a < b< min ( , n/2 + e - a - 1) 

1 U(p + a + l)(a + 1) " V 7 ; 

(36) 

where a = (p/2 + a + l)/(n/2 + e — a) = lim^^oo </>(w), i/ie generalized 
Bayes estimator is minimax. 

Proof. First we prove part (i). Monotonicity of follows from Corollary 13.11 
(ii) since b > min(0, n/2 + e — a — 1). Also using Corollary 13.11 (i) and (1341) 

p/2 + a + l p-2 
< < a = < 2- = c(p, n). 

Hence minimaxity follows from Corollary 12.11 and part (i) follows. 

Next we consider part (ii), the non-monotonic case. The lower bound 
in ([36]) 

-(n + 1)- cM ~ a 



4(p + o + l)(a + l) 
is negative because of ( J35l) . Also Corollary 13.11 (i) and ( l35l) implies 

n v ^ \ p/2 + a + l p-2 

< hm (p(w) =a = — — < 2 — — = c(p,n). 

w^oo n/2 + e — a ra + 2 

Also since the lower bound in fl36l) satisfies 

c(p, n) — a / n/2 + e — a 



- (n + 2) 



4(p + a+l)(a+l) V p/2 + n/2 + e + l 



> _ 1 _f 1 (n + 2)c(p,n; 



a + 1 V. 4(p + a + 1) 
1 ^_ p-2 



a + 1 V p + 2(p/2 + a + l) 

1 2 
a + Ip 
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it follows from Corollary 13.11 (iv) that 



4>{w) < - ^ = Ml 

Additionally, by Lemma I3T41 (iv). 

> (Illiml = _ M2 . 
4>{w) ~ 6+1 

By Corollary 12.21 (ii), it follows that the generalized Bayes estimator is min- 
imax provided 



n + 2) { M 1 -c(p,n)} +4M ^ o 



1 + Mi 

but a straightforward calculation shows that this condition is equivalent to 

b>-{n + 2) Al c M~ a (38) 

Hence the generalized Bayes estimator is minimax since (I36p guarantees (1371) 
is satisfied. This completes the proof. □ 

Note: A recent paper by Wells and Zhou (2008) also derives general- 
ized Bayes estimators in the normal case, some of which have non-monotone 
shrinkage functions <fi (in our notation). The generalized Bayes minimax es- 
timators of Theorem 13.21 (ii) all have non-monotone shrinkage functions and 
their minimaxity cannot be shown by the methods in Wells and Zhou (2008) 
which rely on integration by parts. As noted earlier, the assumption that 
b < causes the usual integration by parts technique to fail. 



4 Generalized Bayes estimators for spheri- 
cally symmetric distributions 

In this section, we consider generalized Bayes minimax estimators for spheri- 
cally symmetric distributions. As shown in Maruyama (2003) and Maruyama 
and Strawderman (2005), the special choice b = —a — 2 in the prior given 
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by leads to the separated joint density of 9 and i], \\9\\ 2 (p/ 2 + a+1 )<r] a 1+e . 
This follows since 



1 / „\ 11/3112 \P/2 / ^ \p/2 

exp 







mil 



v -a-l 



^°V 2+a exp (~\\0\A dt 
oc \\0\\-W+^) v —\ (39) 

if p/2 + a + 1 > 0. Because of this simplification, we make the assump- 
tion b = —a — 2 throughout the rest of this section. Under quadratic loss 
r)(d — 9)'{d — 9), even in the spherically symmetric situation, the generalized 
Bayes estimator is given by E(rjd\X, S) / E(rj\X, S) and hence we have the 
generalized Bayes estimator with respect to our prior, 

Sup J™dv {n+p)/2 ~ a+e f(v{\\X - 9\\ 2 + S})\\9\\~ 2 ^l 2+a+l ^drid9 

Srp J7v {n+p)/2 - a+e f(v{\\x - 0\\ 2 + S})||#||- 2 (p/ 2 + ffl+1 )cV# 

J Rp 9{{\\X - 9\\ 2 + 5)-(n+ P )/2+a-e-l||^||-2( P /2+a+l)^ jj<*> ^(n+^-a+e^ 

f Rp ({||X - #|| 2 + 5)-("+P)/ 2 + a - e - 1 ||0||- 2 (p/ 2+a + 1 )^/ o ^r/("+p)/ 2 - a + e /(r7)^ 

J Rp 9{{\\X - 9\\ 2 + 5)-(n+P)/2+a-e-l||^||-2( P /2 + a+l)^ 

~ f RP ({\\X - 9\\ 2 + S)-( n +p)/ 2 + a - e - 1 \\9\\- 2 W 2 + a + 1 )d9 
if 

/>oo 

rj^+^-^firfdr) < 00. (40) 



Note that this does not depend on / and hence is equal to the generalized 
Bayes estimator in the normal case. In the normal case, as seen in Section 3, 
the estimator is well-defined if a > — p/2 — 1, b > — 1 and e > —p/2 — n/2 — 2. 
Since a = — b — 2, the inequality —p/2 — 1 < a < —1 is also satisfied. 

We note that Theorem 13. 2\ together with the general results of Section 
2, imply that the normal theory generalized Bayes estimators in Theorem 
13.21 remain minimax (but not necessarily generalized Bayes) for the entire 
class of spherically symmetric distributions. The additional restriction that 
b = —a— 2 implies, as noted above, that this subclass is also generalized Bayes 
for the entire class of spherically symmetric distributions. Hence we have the 
following result on the minimaxity of the generalized Bayes estimator with 
respect to ||^|| _2 (p/ 2 + a + 1 )^~ a_1 + e for the general spherically symmetric case. 



15 



Theorem 4.1. (i) Suppose —p/2 — n/2 — 1 < e < —n/4 — 3/2. W7ien 

- P /2-l<a< C(P '" )( " /2+ ; ) ' P/2 - 1 (41) 
7 ~ l + c(p,n) v ; 

where c(p,n) = 2{p — 2)/(n + 2), t/ie generalized Bayes estimator is 
minimax under the moment condition given by (140]) . 

(mJ Suppose e > n/4 — 3/2 and n > 2. When 

-p/2 - 1 < a < a* (42) 

where a* is the larger solution of the equation g(a) = where 

g(a) = (2p + 2n + Ae + 4)a 2 

+ a{2p 2 + 2pn + 12p + 7n + 4(p + 3)e + 10} 
4p 2 + {7/2}pn + 6(p + 2)e + 7ra + 13p + 10, 

the generalized Bayes estimator is minimax under the moment condition 
given by (j4"0l . 



Note: Since a* will be shown to between (—2, —1) in the proof, minimaxity 
of the generalized Bayes estimator with — 2 < a < a* (—1 < —a* — 2 < b < 0), 
which has non-monotone shrinkage factor 0, is new compared to Maruyama 
and Strawderman (2005). 

Proof. First consider that the upper bound of a for minimaxity in Theorem 
, . c(p, n)(n/2 + e) — p/2 — 1 

U*(p,7l,e) = — — r . 

1 + c[p, n) 

A simple calculation gives 

(p-2)(n + 4e + 6) 



u m (p,n,e) + 2 



2(n + 2){l + c(p,n)}' 



Hence u*(p, n, e) is greater than —2 if and only if e > — n/4 — 3/2. To show 
(i), note that when —p/2 — n/2 — 1 < e < —n/4 — 3/2, n, e) < —2 and 
(|4ip satisfies the sufficient condition on a for minimaxity of Theorem 13.21 (i). 
Also since —a — 2 (= b) is nonnegative, the condition on b is also satisfied. 
Hence part (i) follows. 
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Next we show (ii). When e > — n/4— 3/2, the generalized Bayes estimator 
with — p/2 — 1 < a < —2 is minimax because —2 < u*(p, n, e) and —a — 2 (= 
b) > guarantees minimaxity by Theorem 13.21 (i). 

Finally we consider the case where —2 < a < — 1 and e > —n/4— 3/2. The 
condition 6 = —a — 2 < min(n/2 + e — a — 1, 0) in Theorem 13.21 (ii) is clearly 
satisfied under the assumption n > 2. Also a straightforward calculation 
shows that the inequality 

c(p, n) — a . , 4 , 

- n + 2 - \ , -<-a-2 = b 43 

is equivalent to g(a) < 0. Note, (1431) . for a = —1, is not satisfied because 

c(p, n) — a , .c(p,n) — a 

-a-2 + (n + 2 — — —— - -1 + (n + 2 ^ ; - 

v U(p + a + l)(a + l) K J Ap(a + l) 

2(p- 2) 

< -1 + ——7- V 

4p(a + 1) 

< 0. 

Also ( )43|) . for a = — 2, is satisfied because 

civ, n) — a . . civ, n) — a 
-a - 2 + (n + 2 — ; w = in + 2 ; /, > 0. 

Hence g(—2) < and gf(— 1) > 0, which guarantees that the larger solu- 
tion a* of the quadratic equation g(a) = should be between —2 and —1. 
Additionally g(u*(p,n,e)) > because 

civ 77- ) — 

-u»(p, n, e) - 2 + (n + 2)-— — ——7 — — = -u,(p, n, e) - 2 < 0, 

4(p + n, e) + l){a + 1) 

which means that a* is smaller than u*(p, n, e). Hence —2 < a < a* leads to 
minimaxity. This completes the proof of (ii). □ 
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